Recurrences in three-state quantum walks on a plane 
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We analyze the role of dimensionality in the time evolution of discrete time quantum walks 
through the example of the three-state walk on a two-dimensional, triangular lattice. We show that 
the three-state Grover walk does not lead to trapping (localization) or recurrence to the origin, in 
sharp contrast to the Grover walk on the two dimensional square lattice. We determine the power 
law scaling of the probability at the origin with the method of stationary phase. We prove that 
only a special subclass of coin operators can lead to recurrence and there are no coins leading to 
localization. The propagation for the recurrent subclass of coins is quasi one-dimensional. 
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I. INTRODUCTION 

The speed of propagation is affected critically by the 
spatial dimensions in diverse systems, ranging from clas- 
sical statistical phenomena like diffusion or percolation 
[l|, 0] to various quantum transport effects ^] . One of the 
simplest models capturing the essence of classical diffu- 
sive propagation is a random walk on a regular lattice 4] . 
The sensitivity of a random walk to the number of spatial 
dimensions is clearly shown by its recurrence properties. 
A balanced random walker on a regular square lattice re- 
turns to its starting point with certainty in dimensions 1 
and 2 if we wait long enough, whereas there is a nonzero 
probability of escape in dimension 3 or higher, as discov- 
ered by Polya [5| . The probability to return to the origin 
is now called the Polya number. 

The generalization of the discrete time random walk 
to a quantum walk was first introduced by Aharonov, 
Davidovich and Zagury '(|, and independently, from a 
quantum information theoretical perspective, by Meyer 
[7[. In the past few years, there has been a consider- 
able number of studies devoted to quantum walks [8|, 
strongly motivated by their possible application for de- 
signing quantum algorithms [9|. In discrete time walks 
the external, position degrees of freedom are assisted by a 
finite number of internal degrees of freedom, the chirality 
or coin. The situation resembles to the discrete version of 
the Chalker-Coddington model [lfl|. Recent experimen- 
tal demonstration of discrete-time quantum walks with 
many steps in different physical systems (trapped ions 
ll|, optically trapped atoms [I4I, hnear quantum optics 
13j) pave the way to study more complex behavior of 
quantum walks in the laboratory. 

In quantum walks, the question of finding the walker 
involves the measurement of its position. The measure- 
ment unavoidably disturbs the quantum system, which 
has to be taken into account when defining quantities 
to describe its spreading, like the hitting time [ij]. We 
have defined the generalized Polya number [i5| to char- 



acterize the probability of return to the origin of a quan- 
tum walker starting from a localized position. In order 
to minimize the disturbance we have adopted a scheme 
for measurements, where each system from an ensemble 
is observed only once [16|. The quantum Polya num- 
ber exhibits strikingly different behaviour compared to 
its classical counterpart. In dimension 2, the symmet- 
ric quantum walk on a square lattice can be recurrent 
(Polya number 1) or transient (Polya number less than 
1) depending on the coin operator and the initial state 
[17|. For certain coins, localization (trapping) at the ori- 
gin can occur [18[. In dimension 1, although both clas- 
sical and quantum unbiased walks have a Polya number 
1, classically the effect is unstable towards the introduc- 
tion of a bias while in the quantum case there is a region 
of stability fl!!]. These results indicate that the recur- 
rence of walks is not only sensitive to the dimension of 
the underlying lattice, but in the quantum case also the 
coin degree of freedom can ultimately determine its be- 
haviour. In the one dimensional case, one can increase 
the dimension of the coin space from 2 to 3, by allowing 
the walker to stay at its position leading to localization 
\2S\ . Further increasing the dimension of the coin space 
can possibly introduce more interesting effects |2]| . 

The triangular lattice is a planar graph with symme- 
try properties different from the square lattice. It can 
define cither a rank 3 oriented graph or a rank 6 undi- 
rected graph. Quantum walks on triangular lattices have 
been considered for designing effective quantum algo- 
rithms [22]. Discrete-time quantum walks, especially on 
triangular lattices, provide a platform to realize topologi- 
cal phases [231 ■ Propagation of continuous-time quantum 
walks have also been recently studied in [^ . 

In this paper, we consider a two dimensional walk on an 
oriented triangular graph with a three-state coin space. 
The classical symmetric random walk on such a graph is 
recurrent, similar to the two-dimensional square lattice. 
We analyze the quantum walk on this lattice and prove 
that neither localization nor recurrence is possible for 



any unbiased coin, except for trivial cases. We calculate 
the asymptotic decay of the probability at the origin for 
the Grover coin and discuss its dependence on the initial 
state. 

The paper is organized as follows. In Section |lll we 
introduce the three-state quantum walk on the triangu- 
lar lattice and describe the asymptotic behavior of such 
a quantum walk with the help of the method of station- 
ary phase. In Section IIIIl we analyze the recurrence 
properties of the walk driven by the Grover matrix. We 
dedicate Section |lVl to find the general requirements of 
the recurrence in arbitrary coin matrices, we prove that 
only a special subclass of coins leads to recurrence, we 
demonstrate the quasi-onedimensional propagation lead- 
ing to recurrence of such walks by constructing an exam- 
ple. We summarize our results in Section|Vl For the sake 
of completeness we show the recurrence properties of the 
classical walk on the triangular lattice in Appendix [Xj 



II. DESCRIPTION OF THE THREE-STATE 
QUANTUM WALKS 

The Hilbert space of the three-state quantum walk on 
the 2-dimensional lattice has the tensor product structure 



H^Hp(g>Hc 



(1) 



where Hp denotes the so-called position-space spanned 
by the vectors |m) corresponding to the walker being at 
the lattice point m 

•Hp = Span{|m)| meZ^} ^^2^C2). (2) 

The three-dimensional coin-space He has the structure 

He = Span {|e,) I 1 = 1,2,3} , (3) 

where the coin state (chirality) je^) corresponds to the 
displacement of the particle by e^ , a two dimensional vec- 
tor on the plane of walk. The triangle-like topology of 
the walk (see Figure [ij leads to the following form of the 
displacement vectors 



ei = (-^,^), e2 = (l,0), 



63 = (- 



(4) 
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The walks in consideration are unbiased in a sense, that 



E' 



(5) 



A single step of the quantum walk is given by the prop- 
agator 



U = S-{Ip(g)C). 



(6) 



Here Ip denotes the identity operator on "Hp. The oper- 
ator C represents the coin flip and acts only on He- The 
conditional step operator S is responsible for the actual 
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FIG. 1: The lattice points which the walker starting from the 
origin (central vertex) can reach in 3 steps, with the available 
displacement vectors. 



displacement of the walker from position m to m + e^ 
with respect to the coin state je^) and has the following 
form 



5" = ^ ^ |m -I- Gi) (m| ® |e,) (gj 



(7) 



The state of the walker after t time steps is given by 
the application of the time evolution operator © on the 
initial state 

\ij{t)} = ^7/^,(m,i)|m) ® |g,) = C/*|V(0)) . (8) 

The probability of finding the walker at the lattice point 
m at time t is given by the summation over the coin 
states 

3 3 

p{m,t) = ^|(m|®(G,|^(m,i))|2 = ^|V;(m,i)|2 



= IIV-'(m,i 



(9) 



Here, we have introduced the vector of probability am- 
plitudes at the lattice point m 



ip{m,t) = (V'i(m,t),V'2(m,i),'03(m,t))^ 



(10) 



Since we focus on the recurrence properties of the 
three-state quantum walks we consider initial states ac- 
cording to the classical Polya problem, i.e. the walker 
starts localized at the origin. Hence, the probability am- 
plitudes vanishes except at the origin 



V'i(m,0)=0, niT^O. 



(11) 



Nevertheless, we have the freedom to choose the 
initial orientation of the coin which we denote 
by the vector of probability amplitudes tp = 

Since the three-state quantum walk is translationally 
invariant the time evolution equation (|5]) is greatly sim- 
plified by the Fourier transformation 



V;(k,t) = ^V(m,i)e*-"", keK^ 



(12) 



where K = (— 7r,7r]. The time evolution in the Fourier 
picture simphfies into a single difference equation 



V'(k,i) = J7*(k)VXk,0) = U\k)^p. 



(13) 



where ipO*--^ 0) denotes the Fourier transformation of the 
initial state. The propagator in the Fourier picture U is 
given by 



ij = D(k)-C, 



(14) 



where I?(k) is a diagonal matrix determined by the dis- 
placement vectors e^ 

D{k) = Diag (e~**=i , e-'^''^ , e-'^''') . (15) 

The time evolution equation in the Fourier picture (J13p 
is readily solved by diagonalizing the propagator U. Since 
U is an unitary matrix its eigenvalues have the form 

A,(k) = e'''"^"'). (16) 

The solution in the position representation reads 



ipim,t) 



E 



K2 



dk 



g-ik-rrigiwj 



('')*(^;,(k)»^;,(k). 

(17) 
We concentrate on the recurrence nature of the quantum 
walks which is determined by the asymptotic behaviour 
of the probability at the origin. This can be readily an- 
alyzed by means of the method of stationary phase. In- 
deed, the amplitude at the origin reads 






), (18) 



and the probability is simply the absolute square of the 
amplitude. Within the stationary phase approximation 
the important points in the integration domain are those 
where the phase Wj(k) has a vanishing derivative, i.e. 
stationary points. The rate at which the integrals in 
(TTSt decay is determined by the flatness of the phase 
at the stationary points. For a two-dimensional walk 
with a finite number of non-degenerate saddle points the 
probability amplitude at the origin decays as t^^, with 
the probability decaying as t~^ leading to a transient 
walk. A continuum of saddle points (saddle line) leads 
to a probability amplitude decaying at a rate t^^^^, the 
probability decays as t~^ at the origin, thus resulting in 
recurrence. 



III. GROVER WALK 

The Grover operator plays a key role in Grover's search 
algorithm. Used as a coin operator for regular two- 
dimensional quantum walks on a Cartesian lattice it leads 
to the phenomenon of localization for all initial coin 



states except one well defined state. Moreover it is widely 
used in quantum walk based search algorithms. 

The Grover matrix is an orthogonal matrix with ele- 
ments defined as 






(19) 



This matrix has an important symmetry. Indeed, it 
commutes with all permutation matrices. Hence, in the 
d — i case cyclic permutation of the initial chiralities 
will only rotate the probability distribution by 27r/3 in 
the positive or the negative direction. To obtain a ro- 
tationally invariant probability distribution we have to 
choose an initial state which is invariant under cyclic 
permutations. Since the global phase of the quantum 
state is irrelevant we find such a symmetric probability 
distribution results from the initial coin state 



^^5 = ^ (1,1,1)^ 



(20) 



Moreover, the symmetry of the Grover operator implies 
that permuting the initial chiralities does not change the 
recurrence properties of the Grover walk. 

Let us analyze the recurrence of the three-state Grover 
walk in more detail. For that purpose we have to find 
the asymptotic behaviour of the probability at the ori- 
gin. This is determined by the stationary points of the 
eigenenergies Wj(k) of the propagator C/(k) in the Fourier 
picture. For the Grover walk the eigenenergies are given 
by the implicit function 



<D(c^,k) 



det (i/(k) - e*"/) 
sin(fci-|) 



„ / \/3fc2 \ . fki+bj 
2 cos sm 



3sin( 1^1=0. 



(21) 



It is straightforward to show by implicit differentiation 
that the stationary point is ko = 0. Moreover, also the 
second derivatives are vanishing at kg. To clarify this 
statement, we consider the eigenvalues of U{k) for fci = 
or ^2 = —m^i- In both cases we find that one eigenvalue, 
say Ai, is constant and equals unity, i.e. 



UJi{ki = 0,^2) = L0i{ki,k2 = ^=fci) 

v3 



0. 



(22) 



Hence, all derivatives of wi with respect to fci at the line 
^2 = —/^ki vanish, the same applies on the line fci = 
with taking the derivatives with respect to ^2 . Therefore, 
the second derivatives of wi (k) vanish at the stationary 
point ko = 0, i.e. the rank of the Hessian matrix at k 
is zero. In such a case, the method of stationary phase 
indicates that the integrals in the probability amplitude 
(IT51) decay like t~'^/^ as t approaches infinity. Hence, the 



decay of the probability at the origin is given by po (t) ~ 
t-4/3 -^^fi^JQ]^ ig faster than the threshold required for the 
recurrence. We conclude that the three-state Grover walk 
on a plane is transient, i.e. its Polya number is less than 
unity. 

The decay of the probability at the origin can be even 
faster than what we have already found. Indeed, we can 
eliminate the stationary point ko by the proper choice 
of the initial coin state ip. Such ip has to be orthogonal 
to the eigenvector corresponding to Ai at the stationary 
point ko, which is easily found to be 



^'i(ko) 



1 

71 



(1,1,1) 



Hence, the initial coin state has to be of the form 
V'G = (a, ^, c) , a + b + c = . 



(23) 



(24) 



Starting the three-state Grover walk with an initial state 
from the subspace ((24| will lead to fast decay of the prob- 
ability at the origin. 

We illustrate these results in Figures [2] and [S] In Fig- 
ure[2]we choose the initial state tps of (I20|) . The resulting 
probability distribution shown in the upper plot is sym- 
metric and peaked at the origin. However, in contrast to 
the Grover walk on a Cartesian lattice, the current model 
does not exhibit localization. In the lower plot, po(0 is 
shown on a double logarithmic scale. The numerical re- 
sults agree with the power law pQ{t) ~ t~^/^ predicted 
by the method of stationary phase. 

In Figure [31 the initial state was chosen to be V = 



(El 



*47r/3^ which belongs to the subspace 
The peak of the probability distribution at the 
origin vanishes (upper plot). The numerical results of 
the probability at the origin presented in the lower plot 
indicate that the decay rate has doubled to poit) ~ t~^^^. 



IV. RECURRENCE IN THE THREE-STATE 
QUANTUM WALK 

After analyzing a particular case we turn our atten- 
tion to the recurrence behaviour of a general three state 
quantum walk. Let us consider an arbitrary SU(3) coin 
matrix C with matrix elements C^. The characteristic 
polynomial of the propagator in the Fourier picture U(k) 
has the form 

A^* - A2Trf/(k) + A [Ai2(k) + Ai3(k) + A23(k)] -1 = 0, 

(25) 
where we denote 

A,j(k) = e-^''-('=-+^^) (Q,C,, - CuC,,) . (26) 

We have also used the fact that the determinant of f/(k) 
is unity. With the help of the eigenvalues A^ = e*"^ ^^'> of 
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FIG. 2: (Color online) Probability distribution of the three- 
state Grover walk on a plane and the time dependence of 
the probability at the origin (throughout the paper time 
t is measured in dimensionless units, corresponding to one 
step). In the upper plot we display the probability distribu- 
tion after 30 steps. The initial coin state was chosen to be 
ips = -75 (1, 1, 1) . We find that the probability distribution 
is symmetric and peaked at the origin. However, the proba- 
bility at the origin is decaying according to a power law. The 
lower plot is in loglog scale to unravel the scaling law on a 
longer time scale. The straight line corresponds to t~^'^. 



the matrix L'^(k) we can express the characteristic poly- 
nomial (1251) in a different form 



+ A (e^^'^^'-^'>+'^^^^'>^ + gi('-'i(k)+c-'3(k)) _^ gi((^2(k)+cj3(k))A __ 
_ g«(<^i(k)+cJ2(k)+a;3(k)) _ g ^ ^27) 

Comparing the coefficients of the same powers of A in the 
expressions (|25l) and (|27|) we find the relations 



(A°): 



gi(wi(k)+W2(k)+(.J3(k)) ^ I 



(28) 



(A^) : e*("i('')+"2(k)) + gi(wi(k)+t^3(k)) ^ gi(c^2(k)+a;3(k)) ^ 

= Ai2(k) + Ai3(k) + A23(k) (29) 



(A^) 



.i(k) 



aj'^2{k) +gic^3(k) =TrC7(k), 



(30) 



These equations must be satisfied for all k e K^. 

From the relation ([SS)) we could easily replace one of 
the eigenvalues with 



a;3(k) = -a;i(k) - W2(k) . 



(31) 




FIG. 3: (Color online) Probability distribution of the three- 
state Grover walk on a plane and the time dependence 
of the probability at the origin for the initial state tp — 
-j= (1, exp(z27r/3), exp(J47r/3)) . Since t/j belongs to the fam- 
ily (|24|) the behaviour of the Grover walk is changed consid- 
erably. The peak at the origin has vanished, as illustrated 
in the upper plot. Most of the probability is located on the 
ring. Moreover, the decay of the probability at the origin is 
faster when compared to the power law we have found for the 
initial states outside of the subspace (|24p . On the lower plot 
we show the long time behaviour of the probability at the ori- 
gin, thus indicating that the exponent of the power law has 
doubled to —8/3. We use a loglog scale to make the power 
law more visible. The straight line corresponds to t~^'^. 



With this result the equation (15111) simplifies into 



giwi(k) _|_ giw2(k) _|_ g-i(wi(k)+W2(k)) _ Tri7(k) = 

— e Oil + e 022 + e O33 , 



(32) 



where we have expressed the trace of the propagator C/(k) 
explicitly. By multiplying this equation with e' 2 "1 C^) and 
taking the imaginary part we find 



T^i(k) 



IC22I sin ( -k • 62 + Arg(C22) + ^'^iCk) 



sin( -wi(k) 1 = 

|Cii|sin( -k-ei+Arg(Cii) 



Caalsin -k • 63 + Arg(C: 



33 j 



We focus on the recurrence properties of the quantum 
walk, thus the number of the saddle points (ko). Let 
us assume that in a;i(k) we have a continuum of saddle 
points (a "saddle line"). On this line (kg) the gradient 
of wi (k) must vanish in respect of k 



(9k 



0. 



(34) 



k=ko 



One could easily consider that on the saddle line the value 
of a;i(ko) must be constant, thus we could substitute 
a;i(ko) with a constant 20. By taking the gradient in 
respect of k of Eq. ([55| and using the previous assump- 
tions (moving to the saddle line k = ko and substituting 
wi(k) with 2(/) and its gradient with 0) we have 



= 



I + 
I + 



ei|Cii|cos(-ko -ei + Arg(Cii; 
+ e2|C22|cos(-ko •e2+Arg(C22; 
+ eglCasI cos (-kg • ea + Arg (C33) + (j)) . (35) 

The last equation can be satisfied for a continuum of 
saddle points ko (saddle line) only if two of the absolute 
values |Cii|, IC22I, IC33I are zero. It is easy to prove that 
in this case the walk is quasi 1-dimensional. To show 
that let us assume that Cn and C33 are both zero. The 
unitarity of the coin operator introduces two more zero 
elements in the matrix 



Ci3 
C = I C'21 C22 

(^31 C'32 



C" = 



C12 Ci3 
C22 ^23 

C31 



(36) 



(37) 



-0.1 (k) ) . (33) 



The charateristic polynomial of C/(k) with coins C or C" 
only depends on 62, thus the eigenvalues depend only 
on 62. Consequently, the walk propagates only on the 
direction linked to momentum 62 • k, namely ki. In the 
direction k2 the walk decays exponentially, leading to a 
quasi l-dimensional walk. 

It is easy to generalize the method shown above to see 
when the element Cu is not zero and the other elements 
in the diagonal of the coin operator are zero, then the 
walk is quasi 1-dimensional propagating only in the di- 
rection e.i. Note that if we introduce more zeros in the 
off-diagonals of C or C" then the matrices will be permu- 
tation matrices. These permutation coins lead to trivial 
dynamics with no recurrence. 

If the number of the zero elements on the diagonal of 
the coin operator are less than two only isolated saddle 
points. Hence quantum walks with this class of coins are 
transient. 

The last interesting case is when all three diagonal el- 
ements of the coin operator are zero. In this case, the 
equation 1351 has a solution for each kg, i.e. wi(k) is con- 
stant and consequently the walk will be localizing. Nev- 
ertheless, the coin operator is a permutation matrix. We 
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FIG. 4: (Color online) Probability distribution of the walk 
driven by Crec and the time dependence of the probability at 
the origin for the initial state -75(1, 1, 1)^- In the upper plot, 
we show the probability distribution after 30 steps. The cen- 
ter plot shows the quasi l-dimensional behavior of the same 
probability distribution on a density plot with graylevel on 
the log scale. The lower plot shows the t~^ power law scaling 
of the probability at the origin on a loglog scale. The straight 
line corresponds to t~^ . 



encounter here a trivial dynamics consisting in a mere 
relabeling of the position states, in every third step the 
completely relocalized initial state appears at the origin 
independent of the initial coin state. Thus we obtain an 
important result: there are no non-trivial localizing coin 
matrices on a triangular lattice. 

Illustrating the results above we constructed a coin op- 
erator which belongs to the class C 



e(w,k) = dct(u{k)-e"^I 



cos ( ^1 - "^) - V2cos f — 



(39) 



The characteristic polynomial is independent of ^2 , thus 
the eigenenergies are independent of ^2, too. Hence the 
walk is quasi l-dimensional (See Figure |4|). Moreover, 
the double integral in equation p^ simplifies to a single 
integral, leadin g us back to the standard one-dimensional 
quantum walk |l5{ . The decay rate of the probability at 
the origin is t~^ and the walk is recurrent, as seen in 
Figure S) 



V. SUMMARY 

The analysis of three-state quantum walks on a tri- 
angular lattice clearly demonstrates that neither the di- 
mension of the underlying lattice nor the dimension of the 
coin space in itself can determine whether localization or 
recurrence can occur. Although we could construct an ex- 
ample where the walk was recurrent, we have also proven 
that only the coins from the special subclass leading to 
quasi-onedimensional propagation allow for recurrence. 
We have calculated the time exponents of the probabil- 
ity decay at the origin for the Grover walk resulting in 
transience. This behavior is in sharp contrast with the 
recurrence properties of the Grover walk on a regular 
square lattices which traps the walker at the origin with 
finite probability, also known as localization. 

The triangular lattice is one of the few solvable models 
allowing to decide about the recurrence properties of the 
quantum walk. It is certainly interesting to analyze other 
higher dimensional quantum walks from this perspective. 
The presented work is just the fist step in classifying 
quantum walks using the concept of recurrence and a 
number of interesting effects can be expected for related 
models. 
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(38) 



Let us analyze the behavior of a walk driven by Cuec 
observing the eigenenergies given by the characteristic 



Appendix A: The classical 3-way walk on 
2-diniensional Cartesian lattice 



The classical 3-way walk is strongly connected to the 
generalized Pascal's triangle, the so called Pascal's pyra- 
mid [25|, in a similar way as the regular l-dimensional 



classical random walk on a line connects to Pascal's tri- 
angle. The central element in the tth Pascal's pyramid 
is given by 



Coit) 



t\ 



(V3)!3 



(Al) 



This expression is only valid when t is dividable by 3 (as 
the walker could possible return to the origin at every 
3rd step). We normalize the central element, to get the 
probability of the walker returning 



Po(i) ^ ^,C,{t) 



t\ 



3* (i/3)!3 



(A2) 



We approximate the return probability po(0 with Stir- 
ling's approximation 



PoW = 



t\ 



^ 3^3 

3* (t/3)!3 ~ 27ri ■ 



(A3) 



The recurrence properties of the classical random walks 
are represented by the classical Polya number [J] 



7^ = 1- 



S 



where 



5 = ^Po(i) 



(A4) 



(A5) 



In our case po{t) is proportional to i"^, hence the se- 
ries S diverges, the classical Polya number {V) equals 1, 
therefore the classical 3- way random walk is recurrent. 
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